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Фазовая сегментация мультиспектральных 
слабоконтрастных изображений 


Рассмотрены информационные возможности интерференционного метода для сегментации слабоконтрастных 
изображений в фазовом пространстве информативных признаков. Суть подхода основана на формировании 
опорной волны с использованием мультиспектральных составляющих анализируемого изображения, что 
позволяет синтезировать новые фазо-пространственные характеристики и повысить чувствительность анализа. 
Представлены результаты проверки работоспособности метода на модели и изображениях различной физи- 
ческой природы. 


Введение 


Под аномалиями или объектами интереса мы понимаем наблюдения, рассогласо- 
ванные с остальными данными (характеристики отличаются от характеристик соседних 
областей изображения) и обнаружение которых является целью анализа. Сложность их 
выделения, как правило, состоит в том, что, с одной стороны, аномалии во многих слу- 
чаях представляют небольшую область, которую можно принять за шум или дефект 
снимка, с другой стороны, их параметры могут несущественно отличаться от общего 
фона или перекрываться другими объектами. 

Задача выделения визуально неразличимых областей слабоконтрастных изобра- 
жений имеет большое прикладное значение (рентгенограмма, термограмма, биофизические 
снимки и т.д.). В последнее время увеличилось число задач, которые используют ансамбль 
(мультиспектральные/многопараметровые) изображений, например, в таких областях, как 
медицинская диагностика (анализ характеристик ЯМР, томограмм), дистанционное зон- 
дирование земной поверхности (основано на использовании мультиспектральных и 
гиперспектральных методов получения информации), анализ геофизических полей и пр. 

Сегментация подразделяет изображение на составляющие его области или объекты и 
относится к методам высокого уровня обработки [1], являющейся важнейшей стадией 
анализа. Степень детализации зависит от решаемой задачи и конечный успех компью- 
терных процедур анализа изображений во многом определяются точностью сегментации, 
по этой причине значительное внимание уделяется повышению ее надежности. 

Достаточно широко распространены методы анализа на основе использования 
трех спектральных диапазонов с последующим их слиянием в одно результирующее КСВ 
цветное изображение, с последующей сегментацией на основе использования цветовых 
признаков. Однако этот подход весьма редко позволяет идентифицировать слабо- 
контрастные участки из-за незначительного изменения их отражательной способности 
(к сожалению, здесь невозможно привести цветные изображения, демонстрирующие 
этот эффект), поэтому само по себе использование цветовых методов кодирования изобра- 
жений не решает задачи повышения визуальной чувствительности процедуры сег- 
ментации. 
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Ключевой проблемой для осуществления чувствительной сегментации является 
формирование нового пространства признаков на основе анализа локального контраста 
его элементов, который обеспечивают решение поставленной задачи. 

При разработке новых методов анализа изображений возможно использование 
виртуальных аналогов физических методов обработки радиофизических и оптически сиг- 
налов и полей и их математических моделей. 

Целью данной работы является развитие информационных возможностей виртуаль- 
ного интерференционного метода с целью возможности сегментации мультиспектральных 
слабоконтрастных изображений на основе формирования опорной волны с использованием 
мультиспектральных составляющих анализируемого изображения. 


Виртуальное фазовое интерференционное 
преобразование и его характеристики 


В 1935 г. Ф. Цернике разработал метод фазоконтрастной микроскопии, являющийся 
в настоящее время одним из наиболее чувствительных методов анализа слабоконтраст- 
ных оптических изображений. На рис. 1 представлено сопоставление результатов микро- 
биологического изображения на обычном и фазоконтрастном микроскопах (изображения 
были предоставлены доктором И.В. Любошенко, фирма «Р|Базеуле\уу», г. Париж) и схема 
метода виртуальной голографической интерферометрии, описанного в работе [2], ко- 
торый позволяет использовать интерференционные эффекты для повышения качества 
слабоконтрастных изображений. 


а) 6) в) 
Рисунок 1 — Сопоставление результатов обычной (а) и фазоконтрастной (6) 
микроскопии, схема метода виртуальной голографической интерферометрии (в) 


Реализация полного формального виртуального аналога физического метода фазо- 
контрастной микроскопии Цернике для повышения качества и чувствительности анализа 
низкоконтрастных цифровых изображений описана в работе [3] и основана на сле- 
дующих положениях. 

1. Выполнение фазового модуляционного преобразования исходного изображения 
Цх, у) ‚ обеспечивающее переход из пространства действительных яркостей в комплекс- 
ную плоскость и которое предполагает переход от кодирования анализируемых изобра- 
жений в декартовой системе координат к полярной, общий вид которого описывается 
следующим выражением 


ху) 
Ийх : | 
А(х, у) = Ц у)е № = Ве А(х, у)}+ Лас, у} = Аб уе", (0 
где Я - виртуальный аналог длины волны когерентного оптического излучения 
(модуляционный параметр); аргумент Ф(х,у) — угол поворота вектора в комплексной 
плоскости (зависит от соотношения /(х, у) / 4 ). 
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2. Введение виртуального когерентного опорного поля В(х, У) =1, со значением 
и постоянным направлением вдоль оси «х ». 

3. Визуализация синтезированной характеристики суммы этих двух виртуальных 
векторных полей 


[&сь, у) = |А(х,у) + Вбь,у) . (2) 


Виртуальный аналог, в отличие от реального физического метода фазо-контраст- 
ной микроскопии, базирующегося на измерениях только одного параметра — модуля 
суммы векторного сложения, позволяет использовать и другие характеристики [4]. 


В общем случае при |4] *|В| возможно использование четырех характеристик: 
> > 2 > 2 > > > > 2 > 2 > > 
Ко = ПА +В! -2А|В|созФ ; К =.14 +В +24|ВсозФ ; (3) 
АзтФ Ф 
В =аге(К_) = атс +Ф;у=аге(К,)=—; Ф =ехр( 71/1), (4) 
> > 2 
В|-|АсозФ 


ы 


где в в модули, а у, В — фазовые углы векторной суммы и разности 
виртуальных полей соответственно (рис. 2); что позволяет расширить информационную 
базу синтезируемых изображений. 

Кроме этого, возможные модификации метода могут отличаться видом модуля- 
ционного преобразования (1), выбором значения параметра модуляции А и способом 
определения опорного поля. Их определение зависит от характеристик исходного 


изображения, существенно влияет на результат и обусловливает сферу применения. 


А а ЖК 


а) 


Рисунок 2 — Параметры операции векторного сложения (а) и векторной разности (6) 
Как показали результаты наших исследований, с точки зрения задачи сегментации 
слабоконтрастных мультиспектральных/многопараметровых изображений /(х,у)= 


‚= а. целесообразно использование чисто фазового преобразования 


А(х, у) = ехр(/Ф(х, у)) = ехр(/лК(х, у)/ Л). (5) 
Модуль данного векторного преобразования всегда равен единице, а угол его поворота 
Ф(х,у) определяется соотношением яркости конкретной точки изображения Г(х, у) и 


величины А. При выборе значения модуляционного параметра А, равном единице, обес- 
печивается однозначная нелинейная связь между значением яркости /(х, у) пикселя 
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(х, у) и соответствующим ему значением фазового угла Ф(х, у). Данное обстоятельство 


позволяет осуществлять перераспределение соотношения уровней яркостей изобра- 
жения и выявлять слабоконтрастные области, неразличимые исходно. 


При таком подходе в качестве «опорной» волны В\(х, у) возможно использова- 


Е 
ние одного из векторов А, (х, у) ансамбля {1 5х, У) либо сочетание их комбинаций. 

В таком случае для наиболее распространенного трехдиапазонного изображения 
(п=3), которое позволяет использовать цветовое ВОВ кодирование, возможны три 
варианта комбинаций из исходных изображений {1 нЕ 5} для формирования вектор- 


ных полей 
С =П => Г; С =>; С =Б®ЦЬ. (6) 
Каждой паре сопоставимы четыре характеристики : 

В „|= 2 -с0°(Ф, -Ф,)); К, |=/2+с0$(Ф, -Ф,)); (7) 
р, = ат ( В) Л | АФ (8) 

‚ = аг _) = ГС + и) 

. Во - 1- со5(Ф,-Ф,) . 

АФ, =Ф,-Ф,; Ф,, =ехр(/*л*Г,,), 1 =1,2,3;1 = ], (9) 


которые позволяют синтезировать четыре новых цветных ВОВ изображения. 

Варьирование величины А в диапазоне от 0 до 1 соответствует в некотором смысле 
эффекту, подобному изменению величины фокусного расстояния линзы и обеспечивает 
возможность управления визуальной чувствительностью получаемого результата. 
Наилучшая эффективность метода обеспечивается при согласовании значения 4. и ха- 
рактеристик соответствующих областей потенциального интереса (размер и диапазон 
изменения градаций яркостей). 

«Самоорганизующийся» выбор параметра А может быть выполнен на основе 
преобразования с использованием оператора эквализации гистограммы Н исходных 
изображений ( < — стабилизирующий параметр) 


А,(х, у) =ехр(/л [Н (Г, (х, у)) +а]),1 = 12,3. (10) 


Экспериментальные результаты 


Численные эксперименты проводились на ансамблях из трех слабоконтрастных 
изображений различной физической природы (поскольку возможности воспроизведения 
цветных иллюстраций отсутствуют, полученные результаты переведены в палитру Отау). 

На рис. 3 а) представлено модельное КОВ изображение, содержащее семь не- 
различимых человеческим глазом объектов интереса (перепад яркости менее 2%), 
расположенных на плавно изменяющемся фоне. 

Проведенные эксперименты показали, что наиболее информативной, с точки зрения 
задач сегментации ансамбля из трех изображений, является характеристика Д(х, у) 


Ву (х,у) = аг8 (А (х,у)- А(х,у)) ‚1, 1 =1,2,3 1 * ] (0 


и ДФ (9) (разность фаз двух векторов), визуализация которых (рис. 3 6, в) позволила 
выделить все скрытые аномальные участки. Для сравнения на рис. 3 г) представлен 
результат сегментации популярным методом нечетких С-средних на 12 классов с 
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визуализацией по максимуму функции принадлежности, на котором в верхнем левом 
углу с трудом можно различить слабые контуры 3 областей. «Энергетические» хара- 


ктеристики | и | (7) для данной задачи оказались неинформативными. 


(+) 


р. | 
а) 6) в) г) 


Рисунок 3 — Модельное трехдиапазонное изображение: а) оригинал; 
6) В(х, у); в) АФ(х, у) ; г) сегментация методом нечетких С-средних 


На рис. 4 а — в) представлены составные компоненты цветного космического 
изображения окрестности г. Днепропетровска (спутник ЗРОТ). На цветном изображении 
(к сожалению, возможности его привести отсутствуют) превалируют сине-зеленые тона, 
границы отдельных участков достаточно хорошо выделяются, однако структурные осо- 
бенности самих участков являются невидимыми. Рис. 4 г) соответствует изображению, 
синтезированному на основе Д, (х, у) (11), которая позволяет выделить ряд новых, не раз- 


личимых ранее областей (сопоставьте, например, участки в нижней части правого угла 
изображений). Данная характеристика хорошо и, что не менее важно, естественно сег- 
ментирует исходное изображение, о чем свидетельствует связность полученных сегментов. 

Следует отметить, что параметр АФ(х, у) оказался при обработке данных реаль- 


ных изображений менее информативным и в связи с этим далее не приводится. 
Сравнение рис. 4 г), д) и рис. 4 е) наглядно демонстрирует влияние параметра 
модуляции А. Уменьшение значения А (в данном случае с 1 до 0,27) приводит к боль- 
шей детализации — выявлению сегментов с меньшим различием в уровнях яркости 
исходного изображения, а наилучший результат в данном случае был получен с ис- 
пользованием «самоорганизующегося» подхода (10) (рис. 4 е). 


е) 
Рисунок 4 — Космическое изображение: а — в) компоненты КСВ кодирования; 
г) В(х,у) для А = 1,0; д) для д=0,27 
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Интерференционные характеристики являются нелинейными по отношению к яр- 
костям анализируемых ансамблей. Это обстоятельство открывает возможность влияния 
на степень детализации процедуры сегментации за счет инверсии яркостей нормали- 
зованных исходных изображений ансамбля, т.е. возможность использования допол- 
нения ансамбля вида 


Ух, У) =1-Т, „т(хУ), 1=152,3. (11) 
На рис. 5 а) представлено космическое изображение г. Днепропетровска и ха- 
рактеристики В(х,у) (11) для исходного (рис. 5 6) и инвертированного (рис. 5 в) 


изображений, которые свидетельствуют о возможности увеличения чувствительности 
сегментации участков из различных диапазонов яркостей (соответствует реке в данном 


примере). 


Рисунок 5 — Сегментация космического изображения г. Днепропетровска: 
а) оригинал; 6, в) Д(х, у) для исходного и инвертированного изображения, 


соответственно 


На рис. 6 а - в) приведены геофизические поля участка Земной поверхности и ре- 
зультаты их сегментации (рис. 6 г, д) с использованием характеристики Д\(х, у) ‚ для ис- 


ходного ансамбля и его дополнения (11). 


а) 5 в) | 


г) д) 


Рисунок 6 — Геофизические поля: а) гравитационное; 6) электрическое 
потенциальное; в) магнитное; г), д) В(х, у) от исходного и инвертированного 


ансамбля изображений соответственно 
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Данные результаты подтверждают сделанные ранее выводы о высокой информатив- 
ности характеристик В(х,у) для сегментации слабоконтрастных ансамблей из трех 


изображений и о возможности получения дополнительных данных с использованием до- 
полнения исходного ансамбля. На исходных изображениях (включая и результат их КОВ 
кодирования) выделить какие-либо области не представляется возможным. Синтез интер- 
ференционных характеристик Д(х,у) позволяет выделить ряд участков потенциального 


интереса. 


Заключение 


На основании проведенных исследований можно сделать следующие выводы. 

1. Переход от амплитудной к комплексной угловой модуляции яркостей позволяет 
простым и естественным способом осуществлять сегментацию трех диапазонных мульти- 
спектральных изображений с последующим цветовым КСВ кодированием финального 
результата на основе использования фазовых углов векторных разностей синтезирован- 
ных комплексных компонент. 

2. Метод не требует использования какой-либо априорной информации о числе 
и характеристиках областей потенциального интереса. 

3. Варьирование длины волны виртуального когерентного оптического излучения 
обеспечивает возможность тонкой настройки процедуры сегментации, применительно к 
задаче выделения визуально не различимых слабоконтрастных участков. 

4. Возможно обобщение метода на задачу сегментации гиперспектральных изобра- 
жений на основе использования декоррелирующего преобразования ансамбля и первых 
трех «собственных» изображений ансамбля. 
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Г.@. АКйтебйупа, 1.М. Одоуук 

РВазе ЗертещаНоп о? Тьш Мш@зресёга! 512 па|5 

шт Фе ргезеп{ рарег пем даа сарабПиу оЁ фе пуегегепсе 1есбтаче изе4 г зестешайоп оЁ т рВазе 
рога $ ш зрасе оР шогтайопа| ш@1саюг$ \’аз ехаттед. Тве Кеу ра о# 1$ арргоась 15 Базе4 оп е 
теЁегепсе \’ауе огтайоп у Фе аррПсаНоп оР ши Язресвга[ $1па15. ТБаё аПо\уз ю зупезхе пе\м рБазе 
свагасет$Ис$ ап4 ю пиргоуе апа1уз1$ зепууИу. ТВе гези5 оЁ могкше еЁЯсепсу оЁ 1$ пем шефо4 
асшеуе4 Гог се{ат пло4е]$ ап@ 512па[$ аге Четопзгаед. 
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